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In an earlier paper we had developed a method for computing the effective energy transfer between 

any two Fourier modes in fluid or magnetohydrodynamic (MHD) flows. This method is applied to 

a pseudo-spectral, direct numerical simulation (DNS) study of energy transfer in the quasi-steady 

state of 2-D MHD turbulence with large scale kinetic forcing. Two aspects of energy transfer are 

studied: the energy fluxes, and the energy transfer between different wavenumber regions (shells). 

The picture of energy fluxes that emerges is quite complex — there is a forward cascade of magnetic 

energy, an inverse cascade of kinetic energy, a flux of energy from the kinetic to the magnetic field, 

and a reverse flux which transfers the energy back to the kinetic from the magnetic. The energy 

transfer between different wave number shells is also complex — local and nonlocal transfers often 

possess opposing features, i.e., energy transfer between some wave number shells occurs from kinetic 

to magnetic, and between other wave number shells this transfer is reversed. The net transfer of 

energy is from kinetic to magnetic. The results obtained from the flux studies and the shell-to-shell 

energy transfer studies are consistent with each other. 
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I. INTRODUCTION 



In magnetohydrodynamic (MHD) turbulence several scales interact amongst themselves and energy is transferred 
between them. This transfer plays an important role in the generation of magnetic fields. In this paper we study the 
energy transfer in two-dimensional(2-D) MHD turbulence. The MHD and Navier Stokes equations show that energy 
is transferred, in spectral space, to a mode k from modes p and q such that the three wave numbers satisfy the 
condition k + p + q = 0. In an earlier paper (we will refer to this paper as Paper I) |jj we had introduced the idea 
of effective energy transfer between a pair of modes within a triad by the mediation of the third mode, and found the 
formulae for computing these transfers. In this paper we have applied this method to study energy transfer in 2-D 
MHD turbulence. 

In fluid turbulence the dynamics of energy transfer has been well studied. In 3-D fluid turbulence the kinetic 
energy is transferred from large scales to small scales, whereas in two dimensions there is an inverse cascade of kinetic 
energy from small scales to large scales 0. For MHD turbulence there have been various phenomenological 
analytical [6|-[l3|, and numerical studies [fl4|-p3"f to investigate energy spectra, energy cascade rates, etc. Contrary to 
fluid turbulence, the direct numerical simulations (DNS) of MHD turbulence show that the total energy is transferred 
from large scales to small scales both in 2-D as well as 3-D turbulence |lj||24| . There have been theoretical predictions 
of the magnitude and directions of only a few of the various fluxes J^JlJ,^] . 

Magnetohydrodynamic turbulence is sometimes described in terms of the Elsasser variables = u ± b. It has been 
a basic assumption of the phenomenologies for MHD turbulence that the energies associated with z ± are transferred 
to the small scales, and that the flux of the energies is constant in the inertial range in both 2-D and 3-D turbulence 
25j. Thus the the physics of energy cascades appears to be similar in 2-D and 3-D MHD turbulence. 
The magnetic energy of a mode evolves due to two nonlinear terms in the MHD equations — [b.(u.V)b] and 
[b.(b.V)u] — the first exchanges magnetic energy between different scales, and the second exchanges magnetic and 
kinetic energy between different scales. The kinetic energy similarly evolves due to two nonlinear terms — [u.(u.V)u] 
and [u.(b.V)b] — the first one exchanges kinetic energy between different scales and the second exchanges energy 
between the magnetic and the velocity fields. In this paper we have carried out a detailed investigation of the various 
energy transfers in 2-D MHD turbulence. 

Pouquet et al. H studied the energy transfer between large and small scales using 3-D EDQNM closure calculations. 
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They found that the large scales of the magnetic field gain energy in the presence of small-scale residual helicity (which 
is the difference between kinetic helicity and magnetic helicity). The large scale magnetic field, in turn, was found 
to enhance the energy exchange between the small-scale magnetic and the velocity fields. This energy exchange 
resulted in equipartition of small-scale kinetic and magnetic energies. By drawing an analogy between the MHD 
equation for the magnetic field and the vorticity equation for fluids, Batchelor p(ij argued that the transfer between 
kinetic and magnetic energies takes place primarily at small scales. A similar conjecture was made by Pouquet and 
Patterson [^7j . They also conjectured that an inverse cascade of energy from small-scale to large-scale magnetic field 
is the mechanism responsible for the ehancement of large-scale magnetic energy. Contrary to the prediction from 3-D 
EDQNM of Pouquet et al. Q , simulations of decaying turbulence by Pouquet and Patterson |2lj showed that it is 
the magnetic helicity and not the residual helicity which is important for the growth of large-scale magnetic field. 

In a 2-D EDQNM study, Pouquet |(| obtained eddy viscosities for MHD turbulence. She found that the small-scale 
magnetic energy acts like a negative eddy viscosity on the large-scale magnetic energy. The inverse cascade of the 
mean-square magnetic vector potential, and hence the enhancement of large-scale magnetic energy, was conjectured 
to arise due to the destablization of the of large-scale magnetic field by the small-scale magnetic field. She also 
found that the small-scale kinetic energy has the effect of a positive eddy viscosity on the large-scale magnetic energy. 
Ishizawa and Hattori |?| in their EDQNM calculation obtained the eddy viscosity due to each of the nonlinear terms 
in MHD equations. They found that the eddy viscosity due to (b. V)u is positive, leading to a transfer of energy from 
large-scale magnetic field to small-scale magnetic field. In their calculation, energy was also found to be transferred 
from the small-scale velocity field to the large-scale magnetic field. Both, Pouquet 's || and Ishizawa and Hattori's 
calculations give a non-local energy transfer from small-scale velocity field to the large-scale velocity field. 
In a recent work Ishizawa and Hattori [^o| employed the wavelet basis to investigate energy transfer in 2-D MHD 
turbulence. There are some similarities between the observations made by us in this paper and those in the work 
of Ishizawa and Hattori. Using the formalism developed in Paper I we do a thorough investigation of all the energy 
transfers. The wavelet basis used by Ishizawa and Hattori [|o| is a useful tool for investigating spectral properties and 
energy transfer properties in different regions of the flow. In Ishizawa's work, the flow was divided into a turbulent and 
a coherent region which are distinguished by the level of vorticity fluctuations in that region — turbulence is known 
to have high levels of vorticity fluctuations pq] . They found that the energy transfer between velocity /magnetic 
scales occurs more efficiently in the turbulent region and less efficiently in the coherent region. 
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In some of the above mentioned studies |p||27|] a distinction was not made between the energies transferred to 
a mode from the velocity and the magnetic fields Also, the energy transferred into a mode k from different wave 
number regions was not separately computed — only the net energy transfer into a wave number k was obtained. 
The Eddy damped quasi-normal Markovian (EDQNM) closure calculations dealt mainly with coarse-grained energy 
transfer (between large scales and small scales). In another study, Frick and Sokoloff solved a shell model of 
MHD turbulence and calculated only the kinetic energy fluxes between the velocity modes, and the magnetic energy 
fluxes between the magnetic modes. In our simulations we investigate the (1) various energy fluxes arising within 
and between the velocity and magnetic fields, and (2) the fine-grained (considering many wave number shells) energy 
transfer between magnetic and kinetic energies. This gives a more informed picture of the physics of energy transfer 
in MHD turbulence. 

In 3-D MHD turbulent flows, beyond a critical magnetic Reynold number, the magnetic field can grow and reach 
a steady state ]27| , p0| -p2t . This is called the dynamo effect, and is believed to be the mechanism for the generation 
of magnetic fields within astrophysical objects. The non-linear transfer in MHD turbulence is important for under- 
standing the growth process of magnetic energy in a dynamo. Although it is not possible to indefinitely maintain a 
steady state magnetic field in two-dimensional turbulence ]33[ ] , a steady state can be achieved for a finite period of 
time 0. This, coupled with the similarity in physics, allows us to probe the non- linear energy transfer in real MHD 
flows through 2-D simulations. The methods used here are completely generalizable to the 3-D case and can be used 
to study the dynamo problem directly at a later date. 

The paper is organized into the following sections. In Section || we discuss the various energy fluxes and shell-to- 
shell energy transfer rates which were obtained in Paper I and which we will compute in this paper. The numerical 



methods used are outlined in Section [II. In subsection 1IIB, we present a numerical technique to compute the 



fluxes and the shell-to-shell energy transfer rates. The results are presented in Section IV which is divided into three 



parts: Subsection IV A outlines the approach used by us to obtain a steady state; Subsections 1VB and [V C contain 
the observations for the energy fluxes and the shell-to-shell energy transfer rates, respectively. The discussion follows 
m Section [v|. 



II. FORMALISM AND DEFINITIONS 



The MHD equations in real space are written as 
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-^ + (u.V)u = -Vp+(b.V)b + ^V 2 u, (1) 

and 

— + (u.V)b-(b.V)u + Ai V 2 b (2) 

at 

where u and b are the velocity and magnetic fields respectively, and v and /i are the fluid kinematic viscosity and 
magnetic diffusivity respectively. In Fourier space, the kinetic energy and magnetic energy evolution equations for a 
Fourier mode are 

-^-l + 2vk 2 E u (k)= J2 2 5WI(k|p ' q)+ (3) 

' ' k+p+q=0 



]T i^(k|p,q) 



2 

k+p+q=0 



^M + 2 M fe 2 ^(k) = £ ^ M (k| P ,q)+ (4) 

k+p+q=0 



J2 ^(k|p,q) 



2 

k+p+q=0 

where E u (\n) = |u(k)| 2 /2 is the kinetic energy, and E b (k) = |b(k)| 2 /2 is the magnetic energy. The four nonlinear 
terms 5" u (k|p, q), 5" 6 (k|p,q) S* bb (k|p,q) and S^k^q) are 

S""(k|p, q) = -JZe(t[k.u(q)] [u(k).u(p)] + (5) 
i[k.u(p)][u(k).u(q)]) 

S b \k\ Pl q) = -fle(i[k.u(q)] [b(k).b(p)] + (6) 
£[k.u(p)][b(k).b(q)]) 

5" 6 (k|p, q ) = fle(i[k.b(q)][u(k).b(p)] + (7) 
i[k.b(p)][u(k).b(q)]) 

5 bll (k| P ,q) = fle(i[k.b(q)][b(k).u(p)] + (8) 
i[k.b(p)][b(k).u(q)]) 



These terms are conventionally taken to represent the nonlinear transfer from modes p and q to mode k [|2|,|34j . Note 
that the wavenumber triad k, p and q should satisfy the condition k + p + q = 0. The term S'™(k|p, q) represents 



the net transfer of kinetic energy from modes p and q to mode k. Likewise the term S ub (k\p, q) is the net magnetic 
energy transferred from modes p and q to the kinetic energy in mode k, whereas <!5 bu (k|p, q) is the net kinetic energy 
transferred from modes p and q to the magnetic energy in mode k. The term 5 6i> (k|p, q) represents the transfer 
of magnetic energy from modes p and q to mode k. Thus, the quantities S uu (k\p, q), S'" 6 (k|p, q), 5 bb (k|p, q), and 
5 bu (k|p, q) represent the nonlinear energy transfer from the two modes p and q to mode k. We define the following 
quantities : 

r"(k|p|q) = -£e(i[k.u(q)][u(k).u(p)]), (9) 
£ 6,, (k|p|q) = -Re (*[k.u(q)] [b(k).b(p)]) , (10) 
r 6 (k|p|q) = Re (i[k.b(q)][u(k).b(p)]) , (11) 

^"(k|p|q) = Re (i[k.b(q)][b(k).u(p)]) . (12) 

. where Eqs. (p|)-([T^) are the first terms on the right hand side of Eqs. (5)-(8). In Paper I, we showed that the the 
expressions in Eqs. (p[)-(^2[) have the following intepretation : ># u "(k|p|q) is the kinetic energy effectively transferred 
from mode p to mode k; ^ b& (k|p|q) is the magnetic energy effectively transferred from mode p to mode k; ^ tlb (k|p|q) 
is the magnetic energy effectively transferred from mode p to kinetic energy in mode k, and 5$ btl (k|p|q) is the kinetic 
energy effectively transferred from mode p to magnetic energy in mode k. In all these effective energy transfers the 
mode with wavenumber q mediates the transfer but does not itself gain or lose net energy. We discussed in Paper I 
that Eqs. (|9|)-([l^) give us a method of computing energy exchange between any two modes. This method allows us 
to have a much better understanding of cascade rates between scales. 

In this paper, we use the term w-sphere (6-sphere) to represent a sphere in wavenumber space enclosing velocity 
(magnetic) modes. Thus the energy associated with a u-sphere (6-sphere) is kinetic (magnetic) energy. In Paper I, the 
idea of effective mode-to-mode transfer rate was used to define the following energy fluxes and shell-to-shell energy 
transfer rates. The flux 11^5 (K) is the kinetic energy lost by a w-sphere of radius K to the magnetic energy of the 
6-sphere of the same radius, i.e., 

K<(k)= £ J2 ^"( k lp|q); ( 13 ) 

|k|<_ff |p|<K 
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where the summation is restricted such that k + p + q = 0. The flux H£<(K) is the energy lost by a u-sphere to 
modes outside the corresponding &-sphere, i.e., 

n£(iO= E E ^(k|p|q); (14) 

\k\>K\p\<K 

The flux H^y(K) is the energy lost by a 6-sphere to modes outside the w-sphere, i.e., 

Ki(K)= E E ^ 6 (k|p|q). (15) 

|k|>K |p|<Jf 

The flux H%>(K) is the energy lost by modes outside a u-sphere to those outside the corresponding fr-sphere, i.e., 

n£W = E E ^( k lp|q)- (i6) 

|k|>K |p|>K 

These are the only fluxes from u to b modes. The fluxes for the kinetic-to-kinetic and magnetic-to-magnetic energy 
transfers can also be written in terms of the mode-to-mode transfers $ uu and $ bb respectively. Flux H™>(K) is the 
energy lost by w-sphere to modes outside the t/-sphere, i.e., 

K<(K)= J2 E r"(k|p|q). (17) 
|k|>if | P j<if 

Similarly, the flux 11^ (if) is the energy lost by a 6-sphere to modes outside the 6-sphere, i.e., 

n£(*0= E E Ak|p|<0- ( 18 ) 

|k|>K |p|<-R" 

The total flux is defined as the total energy (kinctic+magnetic) lost by the if -sphere to the modes outside, i.e., 

n tot (K) = Ul<(K) + U b <(K) + U^<(K) + U b < > (K). (19) 
Similarly, from Paper I, the shell-to- shell energy transfer rate from the n th it-shell to the m th 6-shell is given by 

T^ = EE^"( k ip|q)- ( 2 °) 

kem pen 

The shell-to-shell energy transfer rate from the n th w-shell to the m th u-shell is given by 

= EE^ U ( k ipi«o- (21) 

kem pen 

The shell-to-shell energy transfer rate from the n th 6-shell to the m th 6-shell is given by 

T n?n = EE^( k IPl f l)- ( 22 ) 
kem pen 



The energy transfer rate from the m 6-shell to the n it-shell, = —T^ n . 

In this paper, we have performed a direct numerical simulation of MHD turbulence and using the values of Fourier 



modes obtained in the simulation, we have calculated various fluxes [Eqs. (|T^)-( 18 )] and energy transfer rates [Eqs. (p0|)- 
(p2|)]. The details of the numerical methods used to compute these quantities and to solve the MHD equations are 
presented in the following section. 



III. SIMULATION DETAILS 
A. Numerical Method 

In our simulations we use the Elsasser variable z ± = u ± b instead of u and b. The MHD equations (1) and (2) 
written in terms of z + and z are 

— T- (B .V) z ± + (z T .V) z ± = -Vp + i/±V 2 z ± + i^V a z =F 

+ (v±/k 2 eq ) V 4 z± + (^/k 2 eq ) V 4 zT + p± (23) 

where — (y ± n)/2. In addition to the viscosity and the magnetic diffusivity in Eqs.(l) and (2), the Eq. fl23| ) 
includes hyperviscosity {v^ jk 2 eq ) to damp out the energy at very high wave numbers. We choose v = fi = 5 x 1CP 6 for 
runs on a grid of size 512 x 512. The parameter k eq is chosen to be 14 in all the runs. The terms F ± are the forcing 
functions. The corresponding forcing functions for the u and the b fields are related to F ± , i.e., F" = (F + + F~)/2 
and F b = (F + — F _ )/2 respectively. In our simulations, we do not force the magnetic field (i.e, F h = 0). Consequently, 
we get F+ = F" = F u = F. 

We use the pseudo-spectral method [[55) to solve the above equations in a periodic box of size 2n x 2ir. In order 
to remove the aliasing errors arising in the pseudo-spectral method a square truncation is performed wherein all the 
modes with \k x \ > iV/3 or \k y \ > iV/3 are set equal to zero. The equations are time advanced using the second order 
Adam-Bashforth scheme for the convective term and the Crank-Nicholson scheme for the viscous terms. The time 
step At used for these runs is 5 x 10~ 4 . All the quantities are non-dimensionalised using the initial total energy and 
length scale of 2tt. For the simulation results shown in this paper a c ~ 0.1. However, we have carried out simulations 
upto a c w 0.9 and the results obtained for higher cr c 's were found to be qualitatively similar to those shown in this 
paper. 
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At each time step we construct F as an uncorrelated random function that is divergence free (i.e., V.F = 0). The 
i-component of F is determined at every time step by 

F x = ( M y/AtFe*, (24) 



, k , 

where T 2 is equal to the average energy input rate per mode, and phase is a uniformly distributed random variable 
between and 2-7T. The y-componcnt of the forcing function is obtained by using the condition of zero divergence, i.e., 

F y = - s/Atfe**, (25) 

We implement the forcing over a k-space annulus 4 < k < 5. The value of T 2 is determined from the average rate 
of the total energy input which is chosen to be equal to 0.1. 

B. Numerical computation of fluxes 

To compute the fluxes we employ a method similar to that used by Domaradzki and Rogallo |j36| . We outline this 
method below using U^(K) as an example. In Eq. ( |l5| ) we substitute the expression for S'" & (k|p|q) [Eq. (|ll| ) ] : 

n%(*)= £ £ tfe(i[k.b(q)][u(k).b(p)]) (26) 

|k|>if \p\<K 

A straightforward summation over k and p involves 0{N 2 ) operations, where TV is the size of the grid, and 
would thus involve a prohibitive computational cost for high N (i.e., high Reynolds number) simulations. Instead, 
the pseudo-spectral method can be used to compute Eq. ( p6| ) in O(NlogN) operations. It involves the following 
procedure. 

We define two 'truncated' variables u > and b < as follows 

!0 if Ikl < K 
(27) 
u(k) if |k| > K 

and 

b( P ) if | P | < k 

h<{ V )={ (28) 

o if |p| > k 



The Eq. (26) written in terms of u > and b< reads as follows 

Kt = Re (*[k.b(k - p)] [u> (k) .b< (p)] ) . (29) 



The above equation may be written as 



U b u <(K)=Re 



5>^>(k)5>(k-p)6<(p) 



(30) 



L k P 

The p summation in the equation above can be recognized as a convolution sum. The right hand side of Eq. ( |30| ) 
can be conveniently and efficiently evaluated by the pseudo-spectral method, using the truncated variables u > and 
h < . This procedure has to be repeated for every value of K for which the flux needs to be computed. The rest of the 
fluxes defined in Eqs. (|l3|)-([l8|) and also the transfer rates defined in Eqs. (po|)-(p2|) are similarly computed. 
In the following section we describe the results of our simulations. 



IV. RESULTS 
A. Generation of Steady State 

The computational time required to obtain a statistically steady state on a grid of size 512 2 is large. So to obtain 
a steady state in simulations on this grid we proceed in stages. We run on a grid of size 64 2 till a steady state is 
achieved. This steady state field is then used as the initial condition to achieve a steady state in a simulation on a 
grid of size 128 2 and so on to grids 256 2 and then 512 2 . 

It is theoretically expected that in 2-D, in the long term, the magnetic energy will decay even if kinetic energy 
is stationary. In our simulations we found that the fields remain steady for the period over which the simulation was 
performed (Fig. 1). Since the N = 512 2 simulations are computationally expensive, we have instead shown the decay 
of magnetic energy for a N — 128 2 simulation. The magnetic energy is found to have a quasi-steady state over the 

time interval extending from approximately t = 20 30 and it decays over the remaining period, as is theoretically 

expected. For the N = 512 2 simulation too, after the period of stationarity, the simulations should show decay of 
magnetic energy. So, strictly speaking, our steady state is only a quasi-steady state. In Fig. 1 we show the magnetic 
and kinetic energy in a quasi-steady state. The Alfven ratio (the ratio of kinetic to magnetic energy) is found to 
fluctuate between the values 0.4 and 0.56. Hence over the steady state magnetic energy dominates over the kinetic 
energy. We compute the fluxes and the shell-to-shell transfer rates over this quasi-steady state once in every unit of 
non-dimensional time. The fluxes and transfer rates are shown in this paper after averaging over 15 time units in 
quasi-steady state. 
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B. Flux studies 



In our numerical simulations we have computed all the energy fluxes defined in Eqs. (|l^)-(jl8|). In this we section 
we discuss these energy fluxes (cascade rates). 

In Fig. 2 we show all the fluxes. The total flux II tot is positive indicating that there is a net loss of energy from 
the if-sphere to modes outside for all K . In the wavenumber region 25 < K < 50, the total flux is seen to be 
approximately constant. This wave number region is the inertial range. 

The net transfer from kinetic to magnetic energy is a sum of the fluxes U^(K), Tl^(K), n^[= —U!^(K)], and 
T¥£y(K). In our 512 2 simulations, the sphere of radius K max — 241 encloses all the modes. We observe from Fig. 
2 that the fluxes ^l^(K max ), H^>(K max ) and H^(K max ) are zero, as there are no modes outside this sphere of 
maximum radius K max . The flux H^(K max ) is found to be positive (see Fig. 2), indicating that there is a net 
transfer from kinetic energy to magnetic energy; since the magnetic modes are not being forced, it is this transfer 
that keeps the magnetic energy constant in quasi-steady state. 

We now describe features of the fluxes II£<, U^< , II*<, II£>, 1T^<, and H b b < observed in our simulations and 
plotted in Fig. 2. We remind the reader that the w-modes within a wave number sphere are called the M-sphere and 
the 6-modes within the wave number sphere are called the 6-sphere. First, we discuss the fluxes that transfer energy 
between the if-modes and the 6-modes. We find that the flux 11^5 (K) is positive — hence, kinetic energy is lost by a 
it-sphere to the corresponding 6-sphere. The flux IE^ (K) is also positive. It means that a M-sphere loses energy to the 
modes outside the 6-sphere. We find that the flux 11^, (K) is negative, which implies that the 6-sphere gains energy 
from modes outside the u-sphere. Thus, all these fluxes result in a transfer of kinetic energy to magnetic energy. 
However, IE^ (K) is negative, implying that there is some feedback of energy from modes outside the 6-sphere to 
modes outside the w-sphere. The net transfer, however, is from kinetic to magnetic. We shall later show that the flux 
(K) plays a crucial role in driving the kinetic-to-kinetic flux (K) . 

The energy gained by the 6-sphcrcs is found to be constant for the approximate range 20 < K < K max . This can 
be seen from the plot of II„< (K) + II^> (which is the negative of the total transfer of kinetic energy to the modes 
within the 6-sphere). The constancy of the flux implies that a 6-sphere of radius K and that of radius of K + AK get 
the same amount of energy from the M-modes. Thus there is no net energy transfer from the u-modes into a 6-shell 
(of thickness AK) beyond approximately K = 20. We therefore conclude that the net energy transfer from u-modes 
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to the &-sphere occurs within the K < 20 sphere. 

We find that there is a kinetic energy gain by the it-sphere from it-modes outside due to the fact that flux n"> (K) 
is negative — this is consistent with the numerical simulations of Ishizawa and Hattori pp| . This behaviour of kinetic 
energy is known as an 'inverse cascade' in literature and is reminiscent of the inverse cascade of kinetic energy in 
2-D fluid turbulence Q and of mean square vector potential in 2-D MHD turbulence (|,|37]]. Note, however that the 
kinetic energy in 2-D MHD turbulence is not an inviscid invariant. Fig. 2 shows that the inverse cascade of kinetic 
energy exists in the wave number range K < 60. Even the modes that are being forced (4 < k < 5) gain energy from 
higher wave numbers u-modes. The source of this energy is the flux Il b u %{K)[= -H%>(K)] which transfers energy 
from higher fe-modes to the higher w-modes and thus effectively forces them. 

We observe from Fig. 2 that there is a loss of magnetic energy from the 6-sphere to the 6-modes outside [see 
IIb> (K)] ,i.e., a forward cascade of magnetic energy — this was also observed in the numerical simulations of Ishizawa 
and Hattori |^0|. Our observation of a forward cascade of magnetic energy contradicts the result of the EDQNM 
closure calculations which predicts an inverse cascade of magnetic energy ||. 

The flux (K) is found to be constant in the inertial range. Using similar reasoning to that given above, we can 
conclude that the net magnetic energy transfer to a wave number shell in the inertial range is zero. However, in this 
case it would imply that the entire energy gained by the shell (of thickness AK) from the ^-sphere of radius K is 
lost completely to the modes outside the region K + AK. Thus the magnetic energy cascades down to higher wave 
numbers independently of the transfers between kinetic and magnetic energy in the inertial range, which are quite 
significant (see below). This flux of magnetic energy is sustained by the fluxes H^(K) and H^(K), both of which 
transfer kinetic energy into the small 6-spheres. 

Figure 3 schematically illustrates the energy fluxes of Fig. 2 for a K sphere of radius K = 20, which is within 
the inertial range. The energy input due to forcing and the small inverse cascade [n^^if)] into the it-sphere from 
higher it-modes, provides the energy input into the it-sphere. This energy is transferred into and outside the 6-sphere 
by H^(K) and H^(K), the latter transfer being the most significant of all transfers (see Figs. 2 and 3). The 
energy transferred into the 6-sphere from the it-sphere [11^5 (If)], and a negligible input from the modes outside the 
it-sphere [— n^S^if)], cascades down to the higher wave number 6-modes. In the higher &-modes, this cascaded energy 
[n^(if)], together with the transfer from the it-sphere, is partly dissipated and partly fed back to the high wave 
number it-modes. This feedback to the kinetic energy is mostly dissipated, though a small inverse cascade takes some 
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energy back into the it-sphere. This is the qualitative picture of the energy transfer in 2-D MHD turbulence. 

The net transfer to each of the four corners of the Fig. 3 sum to zero within the statistical error (which is computed 
from the standard deviation of the sampled data). This is consistent with a quasi-steady-state picture. 

The results presented here for the quasi-steady-state in a forced turbulence remain qualitatively valid even for a 
decaying case — the direction of the various fluxes for the decaying case are identical to that for the forced simulation; 
but as the energy decays, the magnitudes of all the fluxes reduce. 

The fluxes give us information about the overall energy transfer from a wave number sphere or outside-sphere to 
another sphere or outside-sphere. To obtain a more detailed account of the energy transfer, energy exchange between 
the wave number shells are now studied. In the following section we present a discussion on the shell-to-shell energy 
transfer rates in MHD turbulence. 



C. Shell-to-Shell energy transfer-rate studies 

Significant details of energy transfers are revealed by calculating the shell-to-shell energy transfer rates T^, 
T^ n , 7*JJj defined in Eqs. (^l|)-(po|). We partition the k-space into shells at wave numbers k n (n = 1,2,3,...) = 
1, 16, 19.02, 22.62, 2("+ 14 '/ 4 . The first shell extends from k\ = 1 to &2 = 16 — a division of the wave number space 
into smaller shells at the these lower wavenumbers will contain too few modes; the second shell extends from &2 = 16 
to fc3 = 19.02, the m th shell extends from k m to fc( m +i). Thus, the effective shell-to-shell energy transfer rate from 
the n th it-shell to the m th it-shell [Eq. (pl|)] can be written as, 

T r n = E E Er u (k| P |q), (si) 

k m <k<k m+ i fc„<p<fc„ + i q 

and the effective shell-to-shell energy transfer rate from the n th 6-shell to the m th 6-shell [Eq. (|22])j can be written as 

T "nn= E E E^(k| P |q). (32) 

k m <k<k m + i k n <p<k n+1 q 

and the effective shell-to-shell energy transfer rate from the n th it-shell to the m th 6-shell, as defined in Eq. ( |20| ) can 
be written as, 

T ™ = E E E^ U (k| P |q). (33) 

k m <k<k m+ i fc„<p<fe„ + i q 

We categorise the energy transfer between w-shells and the 6-shells as being homologous if the transfers are between 
the corresponding shells (of the same wave number range k n < k < k n+ i, say). Transfers between different shells 
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are therefore non-homologous. Further, transfers (kinetic-to-magnetic, kinetic-to-kinetic, and magnetic-to-magnetic) 
involving shells which are close in wave number space are called local transfers, as is the convention. Two shells are 
considered close if the wave number ratio of the larger shell to the smaller shell is less than 2 (in this study this 
would include the shells between n + 4 to n — 4 from the n th shell). Transfers involving shells more distant are called 
non-local. 

In Fig. 4 we plot the energy transfer rates between it-shells and 6-shells. It is evident from the figure that the 
transfer rates between shells in the inertial range are virtually independent of the individual values of the indices m 
and n, and only dependent on their differences. This means that the transfer rates in the inertial range are self- similar. 
The differences in T^ n for various n are smaller than the standard deviation of the sampled data, indicating that the 
perceived self-similarity is statistically significant. 

We now discuss our simulation results for the non-homologous transfer rates between it-shells and 6-shells. In Fig. 
4 we have shown these transfer rates from the n th it-shell to the m th &-shell by plotting T™, versus m for various 
values of n. We find that for all m, except for m = n — 1 and n, T*^ n is positive. This implies that a it-shell loses 
energy to all the &-shells but gains energy from the (n — l) th and n th ^-shells. The quantity T^ n is found to be small 
for m < n in comparison with the transfer for m > n. Consequently, energy from a it-shell is mainly transferred to 
the ^-shells at higher wave numbers. We had found earlier that the flux H^(K) is much smaller in magnitude to 
the flux H^(K). The fluxes and the transfer rates are consistent with each other. From Fig. 4 we also find that the 
energy transfer rate from the n th it-shell to the I s * &-shell is positive, implying that the I s * 6-shell gains energy 
from the it-shells — this is a nonlocal transfer of energy from the large it-modes to the small &-modes. The flux n„> 



observed in Section IV B (see Fig. 2) is due to this nonlocal transfer. 

The homologous transfers are found to transfer energy from ^-shells to it-shells (see Fig. 4); this is in contrast to 
the net energy transfer and a majority of the non-homologous transfers which are from kinetic to magnetic. We find 
that the energy gained by a it-shell through the homologous transfers is larger than the total loss of energy by the 
it-shell through non-homologous transfers. Consequently, there is a net gain of energy by the it-shells in the inertial 



range. In Section IV B it was shown that the magnetic energy outside a 6-sphere is lost to kinetic energy outside the 
it-sphere. It is now clear that this transfer arises primarily due to the homologous transfers from magnetic to kinetic. 

The Fig. 4 shows energy transfer rates T°^ n , to &-shells from it-shells indexed as n — 5, 6, 7, 8, 9. In Fig. 5 we show 
the transfer rates Tzji from the n — 1 it-shell to all the 6-shells. Recall that the first it-shell comprises the small wave 
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number modes, k = 1 to 16. Comparing the magnitudes of T™f n to the 6-shells from the first w-shell (see Fig. 5) 
and from the w-shells at higher wavenumbers (see Fig. 4), we see that the energy transfer rate from the I s * w-shell 
dominates the transfers from all other u-shells. Hence, there is a large amount of non-local transfer from the first 



w-shell to the 6-shells. In Section [VB we had claimed that there is no net kinetic energy transfer into any inertial 
range 6-shell. Now we show this explicitly in Fig. 6 by plotting the net kinetic energy transferred into a 6-shell 
(=V T^,) versus the shell index m (o in the figure). From the figure, the net energy transferred into the inertial 
range 6-shell can be seen to be nearly zero. We also plot in Fig. 6 the net kinetic energy transfer into a 6-shell from 
all w-shells, except the first (+ in the figure). This quantity is now no longer zero and has a significant magnitude. 
Therefore, the non-local transfer from the first shell plays an important part in balancing the other kinetic energy 
transfers into an inertial range 6-shell. 

We have also computed the transfer rates of magnetic energy T^ n from n th to m th 6-shell. In Fig. 7 we have plotted 
the quantity T^ n versus m for different values of n in the inertial range. We find that the differences in T^ +An ^ n 
for any fixed An(=n — m) are smaller than the standard deviation of their statistical fluctuations computed from the 
sampled data. Hence we can conclude that Tf£ n is self-similar in the inertial range, dependent only on the difference 
An and independent of the location of the shell n. We find that the transfer rates T^ n are negative for m < n and 
they are positive for m > n (see Fig. 7). Hence a 6-shell gains energy from the 6-shells of smaller wave numbers and 
loses energy to the 6-shells of larger wave numbers. Since is self-similar, the energy lost from a shell (n — An) to 
n is equal to the energy lost from n to the shell (n + An) . Thus, the net magnetic energy transfer into any inertial 
range shell is zero, and the energy cascades from the smaller wave numbers to the higher wave numbers. 

We now discuss the simulation results for the kinetic energy transfer rates from the n th to the m th shell. The 
results are shown in Fig. 8 where we have plotted T^ n versus m for various values of n. We find that the most 
dominant transfers are from the n th u-shell to m = n ± 1. From the sign of these transfers we see that kinetic energy 
is gained from (n — l) th shell and lost to (n + l) th shell. This means that the local transfers from the adjacent shells 
result in a forward cascade of energy towards the large wave numbers. The transfers from other shells are largely 
negligible except for the transfer to the shell m = 1 (shown boxed on the left of the figure), which represents a loss 
from high wave number modes to the m = 1 shell. This non-local transfer (to the first shell) produces an inverse 



cascade to the small wave numbers, observed earlier in Section IV B. Now it is clear that this inverse cascade is due 
to non-local transfers. The local and non-local transfers are seen to possess altogether different features: the former 
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is a forward cascade which seems largely self-similar, while the latter is an inverse cascade only to the first shell. 

We schematically illustrate in Fig. 9 the energy transfer between shells. In this figure we show directions of the 
most significant transfers in the inertial range. The arrows indicate the directions of the transfers, and thickness of 
the arrows indicates the approximate relative magnitudes. Since the local transfer rates are self-similar, the transfers 
from any other shells in the inertial range will also show the same pattern. An inertial range 6-shell gains significant 
amount of energy from the smaller w-shells through both local and non-local transfers, and it also locally gains energy 
from the smaller 6-shells (T^J. The energy gained by a 6-shell from the smaller 6-shells is exclusively lost to the larger 
6-shells and the energy gained from the w-shells is mainly lost through homologous transfer (I™J) to the corresponding 
it-shell. A small fraction of the energy is also lost to the larger w-shells. In addition to the energy from the 6-shells, a 
it-shell also gains energy from smaller w-shells by local transfer. The energy of w-shells is mainly lost locally to higher 
6-shells and w-shells, but a significant amount is also transferred to smaller u- modes (T^) by a non-local inverse 
cascade. 

As illustrated in Fig. 9, there is a transfer of energy from the first w-shell to the first 6-shell. This is the most 
significant gain of magnetic energy from the kinetic energy. Under steady state the magnetic energy gained by the 
first shell gets transferred to the higher 6-shells. This transfer yields the forward cascade of magnetic energy discussed 



in Section IV B 

To summarise, we find that there are various types of energy transfers in MHD turbulence: local, nonlocal, ho- 
mologous, and non-homologous. The main local transfers are the forward magnetic energy transfer and the forward 
kinetic energy transfer, the homologous transfer from the magnetic to the kinetic energy (all in the inertial range). 
There is nonlocal transfer from inertial range w-shells to the first 6-shell and the first u-shell, and from the first it-shell 
to the 6-shells. The most significant energy transfer is the transfer of kinetic energy from the first w-shell to magnetic 
energy in the first 6-shell. 

In this section we had extensively described various cascade rates (fluxes) and shell-to-shell energy transfer rates. 
It is clear that the complete picture is quite complex. In fact, some of the features we have seen contradicts earlier 
conjectures and results. These differences are discussed in the next section. 
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V. DISCUSSION 



We have investigated the features of kinetic and magnetic energy transfer between scales at low values of cross- 
helicity in a quasi-steady state of forced 2-D MHD turbulence. Several interesting observations were made in our 
simulations. A summary of the results is given below. For the following discussion refer to Figs. 3 and 9. 

1) There is a net transfer of energy from the kinetic to the magnetic. 

2) There is an energy transfer to the large-scale magnetic field from the large-scale velocity field (w-sphere to 6-shere 
flux 11^5 in Fig. 2), and also from the small-scale velocity field (flux from &-sphere to modes outside the u-sphere 
n^S, in Fig. 2). The former transfer is of a greater magnitude than the latter. The magnetic field enhancement is 
primarily caused by a transfer from the M-sphere to the 6-sphere. Indeed the first few 6-modes get most of this energy. 

3) Other significant transfers, not noted in earlier work, are from the large-scale velocity field to the small-scale 
magnetic field (it-sphere to outside the 6-sphere) and an interesting reverse transfer from the small-scale magnetic 
field to the small-scale velocity field (from modes outside the 6-sphere to modes outside w-sphere). 

4) There is an inverse cascade in the velocity field — this is consistent with the observation of Ishizawa and Hattori 
arising from numerical simulations and is also consistent with the EDQNM closure calculations This 
inverse cascade of kinetic energy is driven by the reverse transfer of energy from magnetic to the velocity field at the 
small scales. Although the flux study points to an inverse cascade of kinetic energy, the shell-to-shell energy transfer 
rates reveal the following feature. There exists both an inverse and a forward transfer of kinetic energy. The inverse 
transfer is primarily nonlocal, coming from the larger wave numbers to the first few modes. The forward cascade is 
local, i.e., between same sized eddies and is self-similar in the inertial range. 

5) There is a forward cascade of magnetic energy towards the small scales. This is consistent with other recent 
numerical simulations pcfl . EDQNM closure calculations also yield a magnetic energy transfer to the small-scales 
0. The magnetic energy transfer is primarily local and is found to be self-similar in the inertial range. 

6) In the inertial range we found a dichotomy between the homologous transfers (defined as the transfers between 
kinetic and magnetic energy shells of same wave numbers) and non-homologous transfers (which are the transfers 
between shells of different wave numbers). These transfers are more complex than had been anticipated |(| and are 
summarised and discussed in detail in Section IV C| (also see Fig. 3). 

The most important outcome of our study can be stated in the following points : (1) Enhancement of magnetic 
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energy occurs due to the energy transfer to the large-scale magnetic field from the large-scale velocity field, (2) there 
is a forward cascade of magnetic energy from the large-scale magnetic field to the small-scale magnetic field. Our 
results are important since they clarify some of the proposals made earlier regarding the physical mechanism for the 
generation of large-scale magnetic field. Pouquet Q had found that the large-scale magnetic energy is destablised by 
the small-scale magnetic energy. Ishizawa and Hattori Q proposed that the large-scale magnetic energy is enhanced 
due to energy transfer from the small-scale velocity field to the large-scale magnetic field. Although we find such 
a transfer, the dominant transfer to the large-scale magnetic energy is from the large-scale kinetic energy. Pouquet 
and Patterson f2?|| proposed in context of 3-D turbulence that there an inverse cascade of energy from small-scale to 
large-scale magnetic field. However, our numerical calculations give a forward cascade in 2-D turbulence. A simulation 
should be done to verify whether this feature is also present in 3-D turbulence. 

The EDQNM closure calculations yield a net transfer of energy to the large-scale magnetic field from the small- 
scales (magnetic -(-kinetic) [BUtI]. In our simulations, the energy lost by the large-scale magnetic field to the small-scale 
magnetic field is much larger than the energy gained by the large-scale magnetic field from the small-scale velocity 
field. Hence, contrary to the predictions of EDQNM calculations, the large-scale magnetic field loses energy to the 
small-scales. A similar observation was also made by Ishizawa and Hattori in their simulations pCj ]. Thus, it appears 
that EDQNM calculations do not yield the correct strengths of the various transfers. 

The picture obtained by us for the energy fluxes and the shell-to-shell energy transfers are consistent and complement 
each other. The fluxes and the transfer rates discussed here could find applications in the dynamo problem. In 
astrophysical objects, like galaxies and the sun, the magnetic field is thought to have arisen due to the amplification 
of a seed magnetic field. Our study has been performed over a quasi-steady state with a low Alfven ratio. In order 
to understand the build-up of magnetic energy starting from a seed value it is important to perform a similar study 
at high Afven ratio. In some of the popular models, like the a-dynamo |^8[ , the mean magnetic field gets amplified 
in presence of helical fluctuations. It must be borne in mind that our caclulations are two-dimsensional and devoid 
of magnetic helicity and kinetic helicity. A three-dimensional calculation (with the inclusion of magnetic and kinetic 
hclicitics) of various fluxes and shell-to-shell energy transfer rates will yield important insights on some unresolved 
issues concerning enhancement of magnetic energy. 

The fluxes also find important applications in various phenomenological studies. For example, Verma et al. ]39| ] 
estimated the turbulent dissipation rates in the solar wind and obtained the temperature variation of the solar wind 
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as a function of solar distance. The various cascade rates discussed here could be useful for various astrophysical 
studies. For example, 11^5 and 11^ can be used for studying the variation of ta of the solar wind. 

The physics of MHD turbulence is still unclear. The studies of various fluxes and transfer rates shed light at various 
aspects which will help us in getting a better understanding of MHD turbulence. 
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FIGURE CAPTIONS 



Fig. 1 Evolution of the total kinetic energy and the magnetic energy for simulations on grid sizes 512 2 and 128 2 . For 
512 2 a quasi-steady magnetic energy is obtained over the period of the simulation. It is demonstrated in the 128 2 
simulation that a quasi-steady magnetic energy eventually decays — a quasi-steady magnetic energy is obtained from 
t = 60 to 100. 

Fig. 2 The simulation results of fluxes. The following fluxes have been plotted : the total flux (II tot ), the kinetic 
energy flux from u-sphere to outside it-sphere njj>, the magnetic energy flux from &-sphere to outside &-sphere 11^, 
the energy flux from it-sphere to &-sphere 11^5 , the energy flux from it-sphere to modes outside the 6-sphere 11^ , the 
energy flux from 6-sphere to modes outside the it-sphere , the energy flux from modes outside w-sphere to modes 
outside the 6-sphere and the net flux out of a 6-sphere (II* < + 11^,) have been plotted in this figure. 

Fig. 3 The schematic illustration of the directions and the magnitudes of the fluxes plotted in Fig. |v| (also see Fig. 3.13 
of Chapter 3) . Also shown are the magnitudes of the kinetic energy input rate due to forcing, and the total dissipation 
of kinetic and magnetic energy. The fluxes are shown for K = 20 but are representative of the entire inertial range. 
All quantities have been time-averaged. The fluctuations of the fluxes (except n^5) and the dissipation rate are 
approximately equal to 0.005. The fluctuations in 11^5 are higher and is approximately 0.01 . 

Fig. 4 The energy transfer rate T^ n from the n th it-shell to the m th 6-shell. The loss of energy from the n th it-shell 
to the m th 6-shell is defined to be positive. 

Fig. 5 The energy transfer rate from the I s * u-shell to the &-shells. 
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Fig. 6 The diamonds (o) represent the net energy transfer into a 6-shell from all the w-shells. The pluses (+) represent 
the net energy transfer into a 6-shell from all zt-shells except the I s * one. 

Fig. 7 The energy transfer rate T^ n from the n th 6-shell to the m th 6-shell. The loss of energy from the n th 6-shell 
to the m th 6-shell is defined to be positive. 

Fig. 8 The energy transfer rate T™ n from the n th u-shell to the m th u-shell. The loss of energy from the n th 6-shell 
to the m th 6-shell is defined to be positive. The boxed points represent energy transfer from the n th w-shell to the 1 st 
zt-shell. 

Fig. 9 A schematic representation of the direction and the magnitude of energy transfer between u-shells and 6-shells. 
The relative magnitudes of the different transfers has been represented by the thickness of the arrows. The non-local 
transfers with the I s * shell have been shown by dashed lines. 
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